Search for i?-parity Violating Super symmetry using Like-Sign 
Dielectrons in pp Collisions at y/s = 1.8 TeV 
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Abstract 



We present a search for like-sign dielectron plus multijet events using 107 
pb~^ of data in pp collisions at y/s = 1.8 TeV collected in 1992-95 by the 
CDF experiment. Finding no events that pass our selection, we set a ■ Br 
limits on two SUSY processes that can produce this experimental signature: 
gluino-gluino or squark-antisquark production with i?-parity violating decays 
of the charm squark or lightest ncutralino via a non-zero \'i2i coupling. We 
compare our results to NLO calculations for gluino and squark production 
cross sections and set lower limits on M{g), M{ii), and M{q). 
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The minimal supersymmetric standard model (MSSM) [|l| is an extension of the standard 
model (SM) that adds a supersymmetric (SUSY) partner for each SM particle and is con- 
structed to conserve baryon number (B) and lepton number (L). The requirement of im- 
parity (Rp) PI conservation is imposed on the couplings: for a particle of spin S, the mul- 
tiplicative quantum number Rp = 1)3-8+^+25 cljgtinguishes SM particles {Rp = +1) from 
SUSY particles {Rp = — 1). If -Rp is conserved, SUSY particles can only be produced in pairs 
and the lightest supersymmetric particle (LSP) is stable. The assumption of Rp conserva- 
tion thus leads to experimental signatures with appreciable missing transverse energy (-^), 
provided that the LSP is electrically neutral and colorless [0. Rp conservation, however, is 
not required by SUSY theories in general and viable Rp violating models can be built 
by adding explicitly B or L violating couplings to the SUSY Lagrangian . Since the LSP 
can be unstable in this case, the standard ^ signature is diluted. 

The results at high from the HERA experiments have sparked interest in ^ SUSY, 
since the excess of events observed at high Q"^ could be explained by the production and 
decay of a single squark: + d ^ q ^ + d, where Rp is violated at both vertices [^j-Q . 
In this scenario, cl (the SUSY partner of the left-handed charm quark) with mass M{cl) — 
200 GeV/c^ is the preferred squark flavor because its associated ^ Yukawa coupling 
less constrained by experiment than the other couplings . Another possibility to explain 
the excess is the production and decay of a first-generation leptoquark; D0 and CDF have 
ruled out this explanation ||TT|| . 

In this Letter, we examine two ^ processes in an MSSM framework that involve the 
same coupling: (1) pp ^ gg (ccl) (ccl) — > c (e^o?) c (e^d) "charm squark analysis"; 
and (2) pp ^ qq ^ ilXi) (qXi) ~^ Q {dce"^) q (dce^) "neutralino analysis". For process (1) we 
assume M(g) > M{g) > M{cl) = 200 GeV/c^, where M(g) denotes the degenerate mass for 
all up- type (except for cl) and all right-handed down-type squarks. The masses of the left- 
handed down-type squarks are calculated using the relations given in Reference 0. These 
assumptions are motivated by the HERA results. Process (2) is a complementary search 
also based on X[2i 7^ 0. It is favored if the size of the coupling is small compared to the 
SM gauge couplings. We separately consider qq production (5 degenerate squark flavors) 
and titi production, and make the mass assumptions: M(xf),M(x2) > M{q) > M(x?), 
where q refers here to either the degenerate squark or ti, and M{xf) ^ 2M(xi). The 
first relation suppresses q xf + ^ ^^'^ the second approximation is generally true for 
most combinations of SUSY parameters, particularly when assumptions leading to gaugino 
mass unification are made. For the case of titi production, we further assume M(xf) > 
M(ti) - M{b) to ensure that Br{ii cx?) = 100% for the relevant case: M{ii) < M{t). 
For these two searches, we make the conservative and simplifying assumption that there is 
only one non-zero ^ coupling. Given the Majorana nature of the gluino and neutralino, 
reactions (1) and (2) each yield like-sign (LS) and opposite-sign (OS) dielectrons with equal 
probability. Since LS dilepton events have the benefit of small SM backgrounds, we search 
for events with LS dielectrons and two or more jets. 

We present results of a search for pp —>■ e^e^+ > 2 jet events using 107 pb"^ of data 
from pp collisions at a center of mass energy of y/s = 1.8 TeV. The data were collected by 
the Collider Detector at Fermilab (CDF) during the 1992-93 and 1994-95 runs of the 
Fermilab Tevatron. At CDF the location of the pp collision event vertex {zyertex) is measured 
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along the beam direction with a time projection chamber. The transverse momenta (pr) of 
charged particles are measured in the pseudorapidity region \ri\ < 1.1 with a drift chamber, 
which is located in a 1.4 T solenoidal magnetic field. Here pr = P sin 6 and 77 = — In tan(^/2), 
where 6 is the polar angle with respect to the proton beam direction. The electromagnetic 
(EM) and hadronic calorimeters are segmented in a projective tower geometry surrounding 
the solenoid and cover the central (|?7| < 1.1) and plug (1.1 < \ri\ < 2.4) regions. A gas 
proportional chamber located at shower maximum in the central EM calorimeter provides 
shower position and profile measurements in both the z and r — (j) directions. 

Dielectron plus multijet candidates are selected from events that pass the central electron 
triggers with Exle) > 9.2 GeV in the 1992-93 run, while for the 1994-95 run there are two 
such triggers, with thresholds of 8 and 16 GeV. The 8 GeV trigger imposes additional 
requirements on the development of the EM shower. In our analysis, we require two electrons 
with Et > 15 GeV. Each electron candidate must exhibit a lateral shower profile consistent 
with that which is expected for electrons, be well matched to a track with pt > Et/2, 
and pass a sliding cut on the ratio of energy in the hadron calorimeter to the energy in the 



EM calorimeter (hadronic energy fraction) At least one electron candidate must also 



pass more stringent identification requirements on its shower profile and hadronic energy 
fraction |TB|. Each electron must pass an isolation cut in which the total calorimeter E^ 



in an 7] — cone of radius R = y(A0)2 + (Ar/)^ = 0.4 around the electron, excluding 
the electron E^, is less than 4 GeV. This helps to remove the background from bb and cc 
production (bb/cc) while retaining much of the sensitivity to the SUSY signal. The ?7 — 



distance ARee = \J{/S.(j)eeY + {ArjeeY between the two electrons must be greater than 0.4 
to avoid shower overlap in the calorimeter. The event \zyertex\ must be less than 60 cm to 
restrict the analysis to the region of the detector that retains the projective nature of the 
calorimeter towers, and both electrons must be consistent with originating from the same 



vertex. Jets are identified in the calorimeter using a fixed cone clustering algorithm ||T6 
with cone size R = 0.7. We require at least two jets with Et > 15 GeV and |?7j| < 2.4, 
separated by ARjj > 0.7, and AR^j > 0.7. Finally, there must be no significant ^ in 
the event: ^ / \JY1, Et < 5 GeV^''^, where Y^Et is the scalar sum of transverse energy in 
the calorimeter for the two electrons and two leading jets. These selection requirements 
are effective in removing the bb/cc and tt backgrounds while retaining the signal. No LS 
candidate events survive this selection, while 165 OS events are retained. 

We calculate the event acceptance using Monte Carlo samples generated with ISAJET 
7.20 |T^, CTEQ3L parton distribution functions |]18|, and passed through the CDF detector 
simulation program. For the charm squark analysis, we examine four values of the gluino 
mass: 210, 250, 300, and 400 GeV/c^ while the charm squark mass, M{cl), is fixed at 
200 GeV/c^. For the neutralino analysis, we create Monte Carlo samples with M{q) in 
the range 100 — 350 GeV/c^. For each M{q), we generate samples for two extremes of the 
neutralino mass: M(x5) = M(q)/2, which corresponds to M{xf) — M(q), and M{xi) = 
M{q) — M{q), the kinematic limit for the decay. 

The dominant SM backgrounds for this search are tt and bb/ cc production, where both 
can give rise to LS ee events. We use ISAJET 7.20 fl^ Monte Carlo samples to estimate 
the sizes of these backgrounds. For ti production and decay, we analyze 25K events (corre- 
sponding to J C dt = 3.3 fb~^) with M(t) = 175 GeV/c^ and a^j = 7.6 pb ||l9l and find zero 
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accepted LS ee events. Top dilepton events typically have appreciable pr and are rejected 
by the ^ significance cut. We study Monte Carlo samples of hb/cc events for two differ- 
ent processes: direct production and final state gluon splitting, and expect a contribution 
of 0.3 ± 0.3 LS events from this source in 107 pb~^. The isolation cut on the electrons is 
efficient in removing this background as semileptonic h quark decays yield poorly isolated 
leptons. The total expected background is therefore consistent with zero events, so we forego 
background subtraction in setting limits. The remaining 165 OS events are consistent with 
the expected contribution of 153.0 ± 14.5 events from SM backgrounds. Drell-Yan produc- 
tion of dielectron pairs accounts for 150.1 ± 14.1 of these events, where we analyze Drell-Yan 
samples generated with pt{Z^ ll*) > 5 GeV/c ||2^ and normalize this production to CDF 
data 1^ before applying the two jet requirement. There is also good agreement in the J2 Et 
distributions for the remaining OS events and for the expected OS background (Figure p. 

The sources of systematic uncertainty on the kinematic acceptances for these analyses 
include initial and final state gluon radiation (ISR and FSR) (4% for the charm squark 
analysis, 4 — 14% for the neutralino analysis), uncertainty on the integrated luminosity (7%), 
electron identification (3%), structure functions (3%), Monte Carlo statistics (1 — 5%), jet 
energy scale (1%), and uncertainty on the trigger efficiency (1%). The total systematic 
uncertainty on the kinematic acceptance is 10% for the charm squark analysis, while for the 
neutralino analysis it ranges from 10% to 16%. 

We set limits on the cross section times branching ratio for the two processes under 
study. In each case we exclude cr ■ Br > Ng5%/{A ■ etrig - J C dt), where iV95% is the Poisson 
95% confidence level (C.L.) upper limit for observing zero events combined with a Gaussian 
distribution for the systematic uncertainty. For both analyses, Nq^% = 3.1 events. The 
acceptance. A, is the product of the kinematic and geometric acceptance and the efficiency 
of identifying two electrons and two jets, and etrig is the trigger efficiency for dielectrons. 
The integrated luminosity is J C dt = 107 ± 7 pb^^. 

For the charm squark analysis, A is a very weak function of M{g) and ranges from 
16.0% to 16.6%. For dielectrons with Erie) > 15 GeV, et„g = 98.4% ± 1.3%. We exclude 
a ■ Br > 0.18 pb independently of M{g). Figure |^ shows the results for the charm squark 
analysis in the gluino-squark mass plane. Exclusion contours at the 95% C.L. are shown 
for two values of the branching ratio Br{cL —>■ ed), where we compare our results to the 
next-to-leading order (NLO) gg production cross section multiplied by the branching 
ratio to LS ee from Reference 0. Our sensitivity vanishes for M{q) ^ 260 GeV/c^. In this 
region is lighter than 200 GeV/c^ (and thus lighter than cl) due to the large top quark 
mass 0, so the decay of ^ bb^ dominates and g ccl is suppressed. Since our analysis 
assumes a non-zero ^ coupling only for c/,, the signal of LS electrons with no ^ disappears 
in this region of parameter space. 

For the neutralino analysis, A is determined for each squark and neutralino mass pair 
and ranges from 3.7% to 15.2%. In this case, etrig = 96.5% ± 1.9%, which is slightly lower 
than for the charm squark analysis because the Et spectrum of the second electron in 
the neutralino analysis is softer. We calculate the upper limit on the cross section times 
branching ratio to LS ee for each squark and neutralino mass combination, and obtain 
a ■ Br limits which range as a function of the squark mass from 0.81 pb to 0.26 pb for 
a light neutralino, and from 0.35 pb to 0.20 pb for a heavy neutralino. Figure |^ shows 
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the results for the neutrahno analysis for the case of titi production. Plotted are our 95% 
CL. upper limits along with the NLO cross section multiplied by the branching ratio 
to LS ee. The branching ratio Br{ti ex?) is taken to be 1.0 p4[. We also assume 



Br{xi — > qq'e) = Br{xi — > ggV) = 1/2, although the actual branching ratios are a function 
of the SUSY parameters ||25|. Since each neutralino decays to e"*" or e~ with equal probability, 
the branching ratio to LS ee is 1/8. The limit is shown for two extremes of the neutralino 
mass, and excludes M{ti) below 120 (135) GeV/e^ for a light (heavy) neutralino. Similarly, 
the results for the case of five degenerate qq production are displayed in Figure^. In this case, 
the NLO cross section [PB| includes a gluino mass dependent t-channel contribution, and we 
assume the branching ratio Br{q —>■ gxi) = l-O- Thus, we set gluino and neutralino mass- 
dependent lower limits on the degenerate squark mass in the range from 200 to 260 GeV/e^. 
The neutralino analysis presented here assumes that the only non-zero ^ coupling is X'i2i- 
Since our analysis does not distinguish the quark flavors in jet reconstruction, however, the 
results are equally valid for any A'^^-;, coupling, for which j is 1 or 2 and /c is 1, 2 or 3. 

We note that our limit for the neutralino decay analysis with 5 degenerate squark flavors 
assumes the branching ratio Br{q qXi) = l-O^ whereas the branching ratio Br{cL ed) 
must be appreciable to explain the HERA results. However, even allowing for Br{q — >■ 
qxT) < 1) our analysis is sensitive to the interesting region of 200 GeV, depending on M{g): 
for example, we can exclude the ^ scenario with Br{q qXi) > 0-^3 for M{g) = 200 GeV. 
For heavier gluino mass, the exclusion becomes weaker. 

In conclusion, we find no evidence for LS dielectron plus multijet events in 1.8 TeVpp col- 
lisions and set a -Br limits on two ^ SUSY processes that could lead to this signature. In the 
charm squark analysis we exclude the scenario of M{cl) = 200 GeV/e^ as a function of M{g) 
and M{q). In the neutralino analysis we set mass limits of M(ti) > 135 GeV/e^ for a heavy 
neutralino (M(x?) = M(ti) — M(c)) and, for the degenerate squark, M{q) > 260 GeV/c^ 
for a heavy neutralino (M(x?) = M{q) - M{q)) and a light gluino (M(^)=200 GeV/e^). 
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FIG. 1. Scalar sum of transverse energy for the two electrons and two leading jets for the 
remaining 165 OS events after all selection (points) and expected background from SM processes 
(histogram). These events are dominated by Drell-Yan e+e" pairs plus two or more jets. 



9 



800 



700 




600 



500 



400 



300 



200 

200 210 220 230 240 250 260 270 280 290 300 

M(g) (GeV/c^) 

FIG. 2. Exclusion region in the g - q mass plane for the charm squark analysis. The branching 
ratio to LS ee is calculated using the scenario in Reference 0], which requires M{g) > M{cl)- 
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FIG. 3. Top: upper limits on the cross section times branching ratio for titi production decay- 
ing to electrons and jets via neutralinos (solid lines). The dashed curve is the theoretical prediction 
for a ■ Br. Bottom: upper limits on the cross section times branching ratio for the production of 5 
degenerate squark flavors decaying to electrons and jets via neutralinos (solid lines). Also shown is 
the theoretical prediction for a ■ Br for three values of the gluino mass: 200 GeV/c^ (dotted line), 
500 GeV/c^ (dot-dashed line), and 1 TeV/c^ (dashed line). 
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